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Abstract 

The effect of bulk dissipation on non critical sandpile models is studied using both multifractal and finite 
size scaling analyses. We show numerically that the local limited (LL) model exhibits a crossover from 
multifractal to self-similar behavior as the control parameters h ext and e turn towards their critical values, 
i.e. h ext — > + and e — > e c . The critical exponents are not universal and exhibit a continuous variation 
with e. On the other hand, the finite size effects for the local unlimited (LU), non local limited (NLL), 
and non local unlimited (NLU) models are well described by the multifractal analysis for all values of 
dissipation rate e. The space-time avalanche structure is studied in order to give a deeper understanding 
of the finite size effects and the origin of the crossover behavior. This result is confirmed by the calculation 
of the susceptibility. 
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1 Introduction 



Sandpiles are prototype out of equilibrium dynamical systems that usually present a self organized 
criticality(SOC). This fundamental concept in modern physics of non-equilibrium phenomena was intro- 
duced by Bak, Tang, and Wiescnfcld [1-5] in order to explain the emergence of scaling behavior and 
fractal structure observed in nature. They proposed the SOC concept as one procedure to describe the 
basic mechanism that creates generic scale-free behavior [1,2]. This behavior emerges when an externally 
driven dissipative system organizes itself into a state where all spatial and temporal events are correlated 
over many orders of magnitude. The main feature of the SOC is that its details are not determined 
by fine-tuning or initial conditions. Moreover, open boundary conditions or bulk dissipation insure a 
balance between input and output flow and allow to non-equilibrium stationary state [6-8]. The SOC 
is postulated be applicable over a wide variety of natural phenomena, spanning from microscopic to the 
astrophysical scale. Bak and his co-workers used the sandpilc [1-4] model as a paradigm because of the 
crude analogy between its dynamical rules and the way sand topples when building a real sand pile. The 
sandpilc model illustrates the SOC for a large class of complex systems [5 and references therein]: biology, 
economics, forest fire models, earthquakes, the game of life, invasion percolation together with systems 
that might be expected to exhibit self-similar and scale invariancc behavior. 

Sandpile systems are modeled as a regular array of columns consisting of cubic sand grains; the usual 
formalization considers each lattice site to be characterized by a state variable h(i), where h is the height 
or local slope of the sand column at a giving site (i). Two mechanisms are crucial in such models: slow 
addition of new grains, which is simply performed by selecting either a random or a fixed column and 
increasing its height by one unit, and the relaxation process. In relaxation process, if the slope in column 
i exceeds a threshold value, then an amount of column sand is redistributed among its neighbors following 
a series of topples which may give rise to an avalanche that subside after a finite period of time. The 
avalanche size 'S' is giving by the number of toppling sites. In SOC, the response to an external pertur- 
bation results in avalanches of all sizes with power-law distributions of the form, D(s, L) ~ s~ T g(s/L D ), 
where L is the system size. The critical exponents r and D, depend on the model one defines. For 
some stochastic sandpile models, the critical exponent r varies continuously [9] by varying the disorder 
in the system. It is believed that the two time-scale separation [6,10,11] (deposition and relaxation) and 
metastability are essential for the existence of scale invariancc in these models. The dynamics of the 
sandpiles have been intensively studied both theoretically [6,12] and experimentally [13, 14], and some 
exact results were derived for abelian sandpile models [15]. Originally, experimental studies showed that 
the sandpile model leads to a clear disagreement with numerical simulations of theoretical models. In- 
deed, if the pile is tilted several degrees above the angle of repose grains start to flow and the system 
exhibits only a first order transition. On the other hand, rice-pile experiments in quasi one-dimensional 
systems display SOC with a power-law distribution of the avalanche size with critical exponent r w 2.02 
[14], and this will occur depending on the detail of the grain level-dissipation. 

Simple rules of the dynamics in sandpile models lead to simple equations according to mean-field 
theory, which is based on the single site approximation of the master equation [6] . One may associate 
rates h ex t and e, respectively, with the addition and dissipation processes. The inverse of the parameter 
h ext is simply the typical waiting time between different avalanches, i.e. w l/h ext . In this case, the 
criticality, which results from non local interactions, is obtained by fine-tuning the control parameters as 
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in continuous phases transitions. 

For a deeper understanding of the significance of SOC a connection to conventional critical points has 
been illustrated through some simple models[16]. As a result, it was shown that the SOC can be under- 
stood as an aspect of multiple absorbing state models since the criticality is reached in the limit h ext — ► 
and e — > + with h ext /e — > similar to Contact Process- like models, and the power-law avalanche dis- 
tribution is found to be a general feature of models with many absorbing configurations [17]. On the 
other hand, using an unified mean field theory the main prediction is that criticality is ensured by the 
divergence of the zero- field susceptibility which is giving by x — dp a /dh ext , where p is the density of 
active sites. In the limit of vanishing control parameters, the stationary state displays scaling that is 
characteristic of SOC. 

The effect of bulk dissipation has been studied in a two dimensional dissipative height sandpile model, 
and it was shown that the SOC behavior doesn't occur for discrete driven dissipative model while it was 
observed for continuous ones given a particular choice of the dissipation rate e[8]. 

Our aim is to study numerically the effect of the bulk dissipation on non critical one dimensional 
sandpile models that obey multifractal analysis, and show that in contrast to the previous dissipative 
systems, the dissipation may influence the avalanche dynamics leading to SOC for LL models. The outline 
of this paper is as follow: in sec. II we define the models, the methods used and explain the problem of 
the multifractality, while sec. Ill is devoted to the study of the dissipation effects on non critical models 
defined in sec. II. In sec IV we give a general conclusion. 

2 Models and methods 

For our systems we assume integer heights h(i) at lattice sites i=l,2...i, where L is the size of the 
system. The local slope of the pile z(i) at site (i) is defined as the height difference between two nearest 
neighbors: 

z(i) = h(i)-h(i + l). (1) 

The boundary conditions of the system are reorganized such that the grains can flow out of the system 
from the right side only. The system consists of a plate of length L, with a wall at i — and an open 
boundary at i = L + 1: 

z(i) = 0, if i = 
h(i) = 0, if i > L. 

The profile of the system evolves through deposition and relaxation. In deposition, a particle is added to 
a random site i and 

h(i) -> h(i) + 1. (2) 

During relaxation, we look at all unstable columns of the pile: a column i of the pile is considered active if 
z(i) > z c . The number of grains N to be toppled is determined by either the limited or the unlimited rule. 
In the case of limited model, the number N is fixed while for the unlimited case, the number of toppled 
grains depends of the slope at the active sites therefore N = z(i)-M, where M is a fixed number. Some 
degrees of non conservation can be introduced in the model by allowing a dissipation of energy during 
relaxation events. In a discrete energy model, one can introduce a probability e that the N transferred 
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particles in the relaxing process are annihilated. The relaxation can take place either according to the 
local (L) or non local (NL) rule 

h(i + 1) = h(i + 1) + AT; (L) 
h(i + j) = h(i + j) + l, j 1....V: {NL) 

Thus four different ID sandpile models [18] may be defined as follow: 

1. Local limited model (LL), 

2. Local unlimited model (LU), 

3. Non local limited model (NLL), 

4. Non local unlimited model (NLU). 

In order to check the effect of dissipation under the total size of avalanches F (which is defined as the 
total number of toppling sites) and the number of grains that drop off the edge, D, we are interested in 
applying two techniques to analyze our data. In the first technique, where the system dynamics exhibit 
self-similarity behavior, the probability distribution function (PDF) p(X, L) of the events F and D is fit 
to the finite-size scaling form represented by, 

p(X,L) = L-Pg(X/L% (3) 

where (3 and v are critical exponents associated with the distribution p and the X quantity, respectively, 
and g is a universal scaling function. This form results from the fractal structure of the system that 
presents self-similar properties. Thus, the fractal analysis is successfully applied to describe the finite 
size effects of such systems which exhibit scale invariance. 

In the other technique, the PDFs are better fit with the multifractal form [18-21] given by the f — a 
representation: 

log w p(X, L)/log 10 (L/L ) = /(Iogio(X/X )/logio(L/io)). (4) 

where Xq and Lq are constants that depend on the value of e and a = log\o{X / Xq) /log\o{L / Lq) . In 
the case where / is a linear function, the probability distribution is then described by only two critical 
exponents. Otherwise, there is a whole spectrum of critical exponents. 

The multifractal analysis is more than a simple fit. The f — a representation is one of the solution of a 
differential equation obtained from a local scale invariance hypothesis [21]. On the other hand, it reflects 
the details of the avalanche structures. Indeed, the study of various space-time structures of avalanches 
in one-dimensional local limited sandpile model [22], without bulk dissipation, shows that there are two 
compact types of space-time avalanches (Fig.l): 

ID linear avalanches, and avalanches with backward events which appear to have a 2D space-time 
structure. The main cause of the lack of self-similarity is the presence of different structures and the effect 
of finite-size, which under ID and 2D avalanches is different. It was found that the study of these two 
types of structures separately gives two different set of critical exponents. As a result, when the scaling 
is studied with the combined effects of both avalanche types, it gives rise to a spectrum of scaling indices 
which reveal the presence of multifractal behavior. The multifractal scaling represented by the equation 
(4) has been found to describe the finite-size effects of the numerically evaluated probability distribution 
function in a better way. 
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Figure 1: The space-time representation of the avalanche for e = 0. The dots represent the active 
sites: (a) example of ID avalanche, and (b) example of 2D avalanche. Consequently, the system exhibits 
a multifractal behavior as it presents two different compact structures. 




Figure 2: The FSS fit of p(D, L) for e = 0.1 and for system sizes ranging from L=60 to L=480. We note 
that the FSS gives a good fit for our data , where vd — 0.02 and (3d — 0.01. 



3 Dissipation effects on non critical sandpile models 

The models described above have been studied for e — [18]. It was shown that they display a 
multifractal behavior. We are interested in the effects of the dissipation on these models, and we study 
whether they can display critical behavior. We conclude by using the two techniques described above 
that for the LL model, the FSS is better for e > e c — 0.04, while the multifractal analysis, represented by 
Eq. (4), gives a satisfactory fit to our data for smaller values of bulk dissipation. However, for the LU, 
NLL and NLU models, only the multifractal fit is appropriate for all values of e. 

3.1 Local limited models 

The local limited models were considered as the simplest ones. During each step of an avalanche, the 
number of grains N which falls to the nearest neighbor is kept constant. Since the critical behavior of 
the system and the avalanche dynamics are independent of the value of N, we will take, for simplicity, 
N = 2. In order to have a good fit, we try to adjust the exponents (3 and v for the FSS analysis, and 
the parameters X and L for the multifractal fit, in such way that for each value of dissipation e, the 
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Figure 3: The FSS fit of p(D, L) for e = 0.1 and for system sizes ranging from L=60 to L=480. We note 
that the FSS gives a good fit for our data , where vd = 0.02 and (3d — 0.01. 

distributions p(F, L) and p(D, L) overlap for different values of the system size L. 

In fig. 2 we present the FSS of the drop number distribution p(D,L) for a dissipation rate e = 0.1. 
However, all curves corresponding to different values of L overlap within the scaling fit (eq.3), with 
(3d = 0.01 and vd = 0.02. In order to emphasize the fact that the FSS gives the best description of 
the finite size effects, we analyse, in fig. 3a, the flip number distribution p(F, L) for different values of the 
system size L. We obtain a good fit for the critical exponents (3p = 0.008 and vp = 0.04. On the other 
hand, the multifractal analysis of p(F, L) (Fig. 3b) shows that the data collapses only for small values of 
F with the exponents F = 1 and L = 1/15, while for F = 25 and L = 1/55 they overlap only for large 
values. Thus, we may conclude that for e = 0.1 the FSS is more adequate than the multifractal analysis. 

In figure 4 the data for both p(F, L) and p(D, L) have been plotted for e = 0.06. We find again, 
as for the previous case which corresponds to e = 0.1, that the FSS form gives an excellent fit over the 
whole range of the data with (3p — 0.21 and vp = 0.26 for the flip number (Fig. 4a), and (3d = 0.01 
and vd = 0.04 for the drop number (Fig. 4c), and the mutlifractal analysis does not give an appreciably 
better fit. It gives a good fit only for intermediate values of the flip number (Fig. 4b) for small values of 
the drop number (Fig.4d) otherwise the data does not collapse for different sizes L. For e = 0.01, fig. 5a 
shows that the FSS analysis works well for small values of F but significant discrepancies arise at larger 
values. There are two different scaling fits which look good over limited range of data. The first one is 
obtained only for small values of F with (3p = vp = 0.38, while the second one is localized for larger 
values with (3p — 1.7 and vp = 1.2. In fig. 5b we have plotted the function log 10/ o(F, L)/log 10 (18L) vs 
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Figure 4: (a)-(c) The FSS and (b)-(d) the multifractal analysis of p(F,L) and p(D,L), respectively, for 
e = 0.06 and for system sizes ranging from L = 128 to L = 512. The FSS gives the best fit of our data 
for all values of e greater than a critical value e c . 

log 10 (_F)/log 10 (18L) for different system sizes L. A good fit is obtained since all the data collapse for 
different values of L. In order to emphasize the fact that the multifractal analysis presents the appropri- 
ate fit, we show respectively in fig. 5c and fig.5d the results of the FSS and the mutifractal analyses for 
the drop number distribution p(D, L). As shown in this figure, the multifractal analysis is a much better 
way to describe the finite size effects than the analysis based on a simple FSS. 

The study of different space-time structures of avalanches for the LL model reveals the existence 
of two classes. For e > e c , the bulk dissipation has a hitchhiker effect on the avalanches. The active sites 
frequently dissipate their energy rather than relax to their nearest neighbor. As a result, the avalanche 
space-time structures are quasi-two-dimensional with very low compactncss(fig.6a-b). Thus, the FSS 
analysis gives the best fit because of the quasi "linearity" of the avalanche structure. For e < e c the 
avalanche structures are more compact than the previous ones(fig.6c-d). The relaxation is more frequent 
than the dissipation thereby leading to more backward events that appears to have a two-dimensional 
structure. Consequently, the avalanche size distribution obeys multifractal analysis. For both cases, 
different structures are generated by varying the dissipation rate e. As a matter of fact, their fractal 
dimension as well as the corresponding critical exponents f3 and v vary continuously. 

From the FSS and the multifractal techniques, we conclude that the LL model with a bulk 
dissipation energy exhibits, at some critical value e c , a crossover from multifractal to self-similar behavior. 
In both regions the critical exponents vary continuously with the dissipation rate e as the space-time 
structure of the avalanches changes with e, producing a low compactness accompanied by a change of the 
finite-size effects. The distribution function of the flip number presents a power law behavior for high 
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Figure 5: (a)-(c) The FSS and (b)-(d) the multifractal analysis of p(F 7 L) and p(D,L) respectively for 
e = 0.01, and for system sizes ranging from L = 128 to L = 512. The best fit is given by the multifractal 
analysis due to the fact that for e < e c , the structure of avalanches became compact and therefore the 
ID and 2D structures begin to react differently. 
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Figure 6: By decreasing the level of the dissipation rate from (a) to (d), the structure of avalanches 
becomes more compact and then the system goes from the situation where the dynamics present a fractal 
structure, to a multifractal behavior. 
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Figure 7: The critical exponent r versus the dissipation rate for e > e c where p(F) presents a power law 
behavior. 




Figure 8: The variation of 2vp + Pf versus the dissipation rate e. The scaling laws have been checked, 
and show that the system presents a non- universal behavior where the critical exponents are e-dependent. 
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Figure 9: The log-log plot of the susceptibility x(e) = dp a /dh ext for a system with periodic boundary 
conditions, and sizes ranging from L = 64 to L = 256. The straight line presents the 1/e behavior. It 
was shown that the susceptibility presents two different variations for low and high values of e. The inset 
delimits the region where the 1/e behavior does not fit with \. 
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Figure 10: The log-log plot of the susceptibility x( e )> f° r different values of the toppling grains N and 
for a system size L = 128. As a result, the dynamics of the system is independent of the details of the 
system. 
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values of e and the critical exponent r increases linearly with increasing the dissipation rate (Fig. 7). The 
FSS technique was used in order to check if there are any scaling laws between different critical exponents. 
As a result, we note that the model presents non-universal behavior since the quantity 2vp + Pf decreases 
linearly with increasing e(Fig.8). A numerical study of the zero field susceptibility X (x = dpa/dh ex t, 
where p a is the density of active sites) of the LL model with a giving value of the toppling grains N, shows 
that it is singular at vanishing values of e, thereby signalling a long-range response function (Fig. 9). The 
critical behavior is recovered in the limit of vanishing driving field, corresponding to the locality breaking 
in the dynamical evolution. For large values of the dissipation ( e > e c ), the zero field susceptibility 
exhibits a power law behavior, 1/e, similar to usual sandpile models. For e — ► 0, it presents a faster 
increase than a simple power law, thereby signalling a more complicated behavior. We mention that for 
e < e c , unlike the usual sandpile models, the density of active sites p a does not vary linearly with the 
external parameter h ext . The behavior of x vs e has been studied for different values of the toppling 
grains number N (Fig. 10), and it's shown that the crossover behavior is insensitive to the details of the 
system in the case of LL models. 

3.2 Local unlimited and non local models 

The investigations of the LU, NLL and NLU models shows that the multifractal analysis is the most 
appropriate technique, for all values of e, and gives a good fit for different system sizes. On the other hand, 
the zero field susceptibility x( e ) does not exhibit 1/e-like power law behavior. It displays rather different 
behavior characterized by some function /(e) that diverges (at vanishing values of e) more quickly than 
any power law behavior. The variation of the number of grains N that flip for any relaxation process, 
and the non locality presented by these models, preserve the non critical power law behavior even for 
high values of the dissipation rate. Since the dynamics of these models under variable and non local 
rules generate back avalanches that are responsible for an appropriate multifractal fit, we believe that the 
dissipation can neither eliminate nor reduce them considerably. Thus, a non power law critical behavior 
of the susceptibility is proposed. The critical exponents that give the appropriate multifractal fit vary 
continuously by varying e. 

4 Conclusion 

To summarize, we have studied the behavior of one dimensional sandpile models with bulk dissipation. 
The various distribution functions of events, for each value of the dissipation rate e, depend on the system 
size. However, an investigation of the scaling properties of the one dimensional LL, LU, NLL and NLU 
models have been established using both simple FSS and multifractal analyses of our data. Using both 
FSS and multifractal techniques, we have shown that for the LL model the first one works better for high 
values of the dissipation rate e > e c , while the second gives a more appropriate fit for low values of e. The 
analysis of the space-time structure of avalanches shows that such structures are quasi-two-dimensional 
for high values of e, whereas they become more compact 2D areas for low values. Consequently, the self- 
similarity breaks down for e < e c , leading to multifractal behavior resulting from two kinds of space-time 
structures (ID and 2D compact avalanches) which start to react differently. In order to emphasize the 
result obtained using the techniques mentioned above, the crossover behavior has been checked by the 
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calculation of the zero field susceptibility x(e) in view of the fact that it exhibits 1/e -like power law 
behavior for large values of e, while for small values such behavior is replaced by a rather different one. 
We note that since the crossover behavior is localized at very low values of e, and that the system deviates 
slowly from the power law behavior, the susceptibility, the FSS and multifractal techniques do not allow 
high precision of the estimated value e c = 0.04. 

The LU, NLL and NLU models do not display any crossover behavior, and they preserve the multifractal 
analysis since their intrinsic dynamics produce frequent back avalanches (for all values of e) that are 
responsible for such behavior. The critical exponents obtained either within the multifractal analysis or 
the FSS vary continuously with the value of the dissipation rate e. 
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